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Abstract 

Consider the motion of a material point of unit mass in a central field determined 
by a homogeneous potential of the form (— 1/r"), a > 0, where r being the distance 
to the centre of the field. Due to the singularity at r = 0, in computer-based simu- 
lations, usually, the potential is replaced by a similar potential that is smooth, or 
at least continuous. 

In this paper, we compare the global fiows given by the smoothed and non- 
smoothed potentials. It is shown that the two flows arc topologically equivalent for 
a < 2, while for a > 2, smoothing introduces fake orbits. Further, we argue that 
for a > 2, smoothing should be applied to the amended potential c/(2r^) — 1/r", 
where c denotes the angular momentum constant. 

Keywords: central field, singular homogeneous potential, smoothing, regularized 
vector field, topological equivalence 



1 Introduction 

For large particle systems, a principal tool of investigation is computer- 
based simulation. In a variety of problems the interaction of the particles is 
determined by a potential that is undefined at collisions. A common technique 
in dealing with the vector field singularities is to replace the potential with a 
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smooth, or at least continuous, function. This procedure is called smoothing, 
or, in physics terminology, softening. 

Smoothing was introduced in 1963 by S.J Aarseth cf. [1], [2], in the context 
of numerical simulations of galaxies. Since then, smoothing has became a 
commonly used technique in numerical modeling of large particle systems (see 
for instance, [6], [10], [9] or [11] ). 

Understanding the modifications induced by smoothing in large particle sys- 
tems still remains a challenging task. A first step is to look at systems formed 
by two particles, but even in this simplified context, one is faced with difficul- 
ties; see, for instance, the analysis presented in [5], where several conjectures 
concerning the convergence of the approximation methods are stated. 

Closely related to smoothing is the concept of regularization: they both tar- 
get singularities in the fiow as induced by the singularities in the vector field, 
but the resolution is different. Smoothing modifies the vector field. Regular- 
ization relies on a qualitative analysis of the phenomena near singularity and 
is achieved in two distinct steps. First, new parametrizations are applied, both 
time-dependent and -independent, leading to a regularized vector field, that 
is a vector field free of singularities. The phase space in the new coordinates 
is extended to include the singularity set, now blown-up into a physically fic- 
titious and invariant manifold, usually called the collision manifold. Second, 
analysis of the fiow on the extended phase space is performed in order to 
decide whether solutions asymptotically reaching the collision manifold can 
be matched to solutions asymptotically leaving the collision manifold, while 
preserving good behavior with respect to initial data. If such a matching is 
possible, then the fiow may be extended (at least continuously) to include 
orbits ending/ starting in collision. When this extension is performed, then the 
problem is said to be regularized. (For more on regularization, see [8] or, from 
a more physical point of view, sec [12].) 

We also mention the paper of BcUcnttini et al. [4], where regularization 
is seen from the different perspective of approximating collision solutions by 
solutions of the smoothed fiow. While analyzing a system where the interaction 
is given by homogeneous potentials of the form — 1/r", a > 0, the authors 
convey that their procedure leads to a larger set of regularizable problems 
than in the standard treatment. Moreover, the smoothing chosen is irrelevant, 
as long as it provides a fiow free of singularities. 

In this paper we question the appropriateness of smoothing when motion 
both near and far from collision is under scrutiny. Our analysis is performed 
within the class of homogeneous potentials to which a standard potential 
smoothing 

t/,(r;a):=-l/(r2 + e2)"/2, e > 0, (1) 

is applied. Within negative energy levels, we focus on the topological equiv- 
alence of the non-smoothed and the smoothed fiows outside the collision set. 
We show that for a > 2 the two flows are not topologically equivalent and 
thus smoothing of the form (1) generates orbits that do not correspond to 
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orbits of the real non-smoothed motion. For this case, we introduce the idea 
of smoothing the amended potential and show that, with such a modification, 
the two flows are topologically equivalent. 

Employing a technique similar to that of McGehee [8] , we choose to describe 
the dynamics in a parametrization where the non-smoothed flow is nonsingular 
and where the phase space is extended to include the collision set, now blown- 
up into a one-dimensional manifold (this is the first step of regularization 
as described above). The orbits lie on compact three dimensional manifolds 
which are level sets of the energy integral for negative energy values. Since 
the regularized vector field preserves the SO{2) equivariance of the original 
problem, dynamics can be studied in a reduced three dimensional space. Here 
orbits can be easily visualized as curves determined by the intersection of 
the two surface integrals, the energy and the angular momentum. This allows 
us to compare the orbital pictures of the non-smoothed versus the smoothed 
problem. 

While the non-smoothed reparametrized fiow includes the orbits on the 
collision manifold, our analysis refers only to orbits outside of it and concerns 
only orbit topology. We do not refer to regularization of solutions (the second 
step of regularization, as outlined above) and we do not focus on issues related 
to approximating the non-smoothed solutions by smoothed ones. 

The paper is organized as follows: we begin by briefiy reviewing known 
facts about dynamics of two particle systems. Next, we reparametrize the vec- 
tor fields of the non-smoothed and smoothed problems such that the collision 
set of the non-smoothed problem is blown-up into the aforementioned collision 
manifold. Using the SO{2) symmetry to reduce the phase space to three di- 
mensions, we study relative equilibria and examine symmetries of the reduced 
fiow. Further, we analyse and compare of the orbits of the non-smoothed and 
smoothed fiows, drawing the conclusion that for a > 2 the two fiows are not 
topologically equivalent. In the last section we argue that for a > 2, smooth- 
ing should be applied to both the potential and the rotational non-inertial 
term, leading to the idea of a smoothed amended potential. Moreover, we 
show that when such a smoothing is applied, the topological equivalence of 
the non-smoothed and smoothed fiows outside the collision set is achieved. 



2 Equations of Motion 



Consider the two degree of freedom Hamiltonian system given by the system 
of first order ordinary differential equations: 



where q = (x, y) e R^, p = {Px,Py) £ R^- The function 

^e(q,p) :=y + C/e(q;a) 
is the Hamiltonian of the system and is a "smoothed" potential 

where e > is a parameter and r = |q| = \/x^ ^ y^. For e = the potential 
t/o(q;Q;) reduces to the classical homogeneous potential, in which case the 
vector field defined by (2) has a singularity at q = 0. 

Since the system is Hamiltonian, it is well known that the total energy is 
conserved. Consequently, the level sets of if^ are invariant under the fiow of 
(2) 

-f^e(q(^))P(^)) = h = const. 
Due to radial symmetry, the angular momentum is conserved as well and we 
have: 

x{t)py{t) - y{t)px{t) = c = const. 
Therefore the system has two independent first integrals in involution and it 
is integrable by the Liouville-Arnold theorem. 

Since t/e(q, a) : ]R^\{0} — M is real analytic, standard results of differential 
equation theory guarantee, for any initial data (qo,Po) = (q(0),p(0)) G \ 
{0} — > M, the existence and uniqueness of an analytic solution defined on a 
maximal interval [0,t*), where < < oo. If < oo, the solution is said to 
experience a singularity. 

For e = 0, the singularity in the vector field induces singularities in the 
solution. Singularities specific to particle systems are given by collisions, which 
occur when q(t) — > as i — > i*. In [8], McGehee showed that coUisions are 
the only possible singularities of (2). He also proved that if a e (0,2), then 
(qo, Po) leads to a collision if and only if the angular momentum is zero (that 
is qo X Po = 0) and that if a > 2, then the set of initial conditions leading to 
a collision is rather large, including for ck > 2 zones where the energy integral 
is negative. 

The standard methodology in dealing with collisions in n-particle simula- 
tions is smoothing the potential by setting e > 0. Then the vector field (2) is 
real analytic for all initial data (qo, po) G x M^, and the associated system 
of differential equations admits a unique global analytic solution. 

3 The flow of the smoothed potential 

3. 1 Topological description of the energy surfaces 

From now on, unless otherwise stated, the energy h is assumed to take 
negative values. 
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Using a technique similar to McGehee [8] we consider the following trans- 
formations {r,Pr,0,pg) {r,v,9,u) defined by 



a-2 

u = pgr 2 . 

V — PrV^ . 



(3) 



This transformation is a diffeomorphism from (0, oo) x R x x R to itself 
(where is the unit circle). Further, we rescale the time parametrization by 



a+2 

dr — r 2 dt. 



(4) 



It is useful to keep in mind v and u are re-parametrised linear and angular 
momenta, respectively. The equations of motion take the form 



r — rv 
V 



„.2 , a.,2 _ ar"+2 



(5) 



U 



'a-2- 



where prime denotes differentiation with respect to the independent variable 
T. In the new coordinates the conservation of energy integral reads 



2 I 2 
U +V — 



2r' 



+ 2r"\h\ = 0. 



(6) 



(r2 + e2)"/2 

Prom above we deduce that the phase space is foliated by the energy surfaces: 
T.h,e := {(r, V, ^, li) e [0, oo) X R X §^ X R I F,(r, v, 6, u) = 0}, 

where 



F^{r,v,9,u) := u'^ + v'^ - 



2r° 



'j.2 _|_ £2^a/2 



+ 2r"|/i|. 



For e = 0, the transformations (3) and (4) have important consequences. 
First, the dynamics given by (5) has no singularity at r = 0; in fact in the 
new coordinates the system extends analytically to all of (r, v, 9, u) e [0, oo) x 
R X X R, and hence the equations of motion are regularized. Second, the 
submanifold r = is now invariant under the fiow. Third, the energy relation 
also extends to the set r = giving 



+ -y^ - 2 = 0. 



Let 



Mo :={(r,'f;,^,ii) |r = 0, ^ e §\ + - 2 = 0}, 



(7) 
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that is Mo is the boundary of the extended phase space. Note that, since each 
energy surface has the same boundary r = 0, the set A/q is independent of 
h. Topologically, it is a two-dimensional manifold embedded in x and 
it is diffeomorphic to a torus. We call Mq the collision- ejection manifold or, 
simply, the collision manifold. 

In the new paramctrization, orbits which previously reached r = in a 
finite time now tend asymptotically toward collision manifold. Further, orbits 
which previously passed close to collision now spend a long time near Mq. By 
continuity with respect to initial data, the flow on collision manifold, although 
lacking physical meaning, provides useful information about collision and near- 
collision solutions. 



For e > 0, using the energy relation (6), we deduce that on each E/^ g the 
radial coordinate is bounded by: 



< r < i?^ 



n\h\. 



2 1 a 

I - (8) 



In physical space this means that the maximal distance in between points is 
bounded by the value i?™"^. 

Remark 3.1 Inequality (8) is meaningful provided e is small enough, that is, 
Henceforth we will assume this condition is always satisfied. 



The energy surface T,h,o is a smooth manifold. For e > 0, if a < 1, 'Eh,e is 
also a smooth manifold; if a > 1, the energy surface is not smooth at (0, 0, 9, 0) 
since VFg at (0, 0, 9, 0) = (0, 0, 0, 0) as it can be readily seen from 



dF, / -2a 2r^a 



dr 



Proposition 3.2 For e > and negative energy T,h,e is diffeomorphic to x 
S-*^ for a <1 and homeomorphic to it for a> 1. 

Proof In the space (r, u, v) e [0, oo) x M^, H^{r, v, 9,u) = describes a surface 
of revolution. This surface is orientable, compact, connected and of genus 0. It 
is smooth for a < 1, and thus it is diffeomorphic to a two-sphere S^. It is not 
smooth for a > 1 and therefore it is homeomorphic to a sphere. Consequently, 
for a < 1 (a > 1), E/j^g is diffeomorphic (homeomorphic) to x S^. □ 

In the new coordinates, the conservation of the angular momentum trans- 
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lates into the presence of invariant surfaces of the form 

^ _ i{{r,v,e,u) e[0,oo) xRxS^ xR\ur^ = c} if a<2 
\{{r,v, 9, u) e[Q,oo) xRxE>^ xR\u = cr"^} if a > 2. 

Taking into account energy conservations, we deduce that the phase space is 
fohated by invariants sets obtained as intersections 12h,e H Fc. 

3.2 Reduced dynamics 

We now return to the system of equations (5). In order to study the dy- 
namics it is convenient to exploit the fact that 6 does not appear exphcitly 
in the equations. This allows us to reduce the four-dimensional phase space 
by factoring out the flow by S^. The three dimensional reduced phase space is 
described by (r, v, u) G [0, oo) x R x R and a vector field given by: 

r' — rv 

-'-^' + - (^^S^ (10) 



Defining 
the energy relation (6) reads: 



u' = {'^)uv. 



u' + v' = Ur). (11) 



With a slight abuse in notation, we choose to call the energy surfaces in the 
reduced space 

E^,, := {(r, v,u) e[0,oo) xRxR\u^ + - fe{r) = 0}. (12) 
Likewise, the angular momentum invariant sets given by (9) are: 

^ _ U{r,v,u) e [0,oo) xRxR\ur^ ^ c} if a<2 
~ [{(r,^,^) e [0,oo) X R X R I K = cr"^} if a > 2. 

3.3 Relative equilibria 

We now study the existence and nature of the equilibria of the reduced 
dynamics. Such fixed points which, in addition, are located outside the collision 
manifold (that is, those with r > 0) correspond in the unreduced space to 
relative equilibria. In physical space, they represent uniform circular motions 
around the center of mass. 
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Case e = 0. Fixed points with r = are found on the coUision manifold Mq 
and are situated at C± = {0,±\/2,0). Direct calculations show that they are 
saddles. In physical space they correspond to radial fall to or ejection from 
Mq. The radial coordinate of fixed points with r > is given by: 

2\h\r'^ ^{2-a). 

This equation has a unique solution 

for a < 2 and no positive solutions if a > 2. Thus relative equilibria are 
present only for o; < 2 and they are of the from Rq — (rg, 0, iug), where u^, 
may be determined by substituting r = Tg and v — into the energy relation 
(11). A direct verification shows that they are centres. 

This is not surprising, since it is well-known that motion in an attractive 
potential of the form 1/r", a < 2 possess such equilibria in the reduced space 
(corresponding to circular orbits in the physical space) as long as the angular 
momentum is non-zero (see [7]). 

Case e > It is immediate that the origin O := (0, 0, 0) is a fixed point. The 
remaining solutions, if they exist, are of the form Rf — {r^,0,±u^). Setting 
V — in the second equation of (10), we obtain 

j.a+2 

- °^(^2 + ,2)f+l = 

and, since by the energy relation 

- Mr) =u'- + 2r"|/i| = 0, 

-I- e^j 2 

we have to solve 

j,a+2 2r" 

2^ = ~ ( 2_L + 2r"|/i|. 

{r^ + e''')'^^ (r"' H- e^j 2 

This amounts to studying the zeroes of the function r"" jir^ -|- e^)^'^/^+^)/i(r) 
where 

h{r) ■= (2 - ay + 2e2 - 1\h\y + e2)f +\ 

If e is sufficiently small, /i(0) = 2e^(l — j/'le") > and hir) — > — oo as r — oo. 
This implies that h{r) has at least one zero for r > 0. The derivative is 
/i'(r) = 2rg{r) with 

g{r) := (2 - a) - 2|/.|(| + l)(r^ + 6^)t. 
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Fig. 1. The reduced orbit space for e = and a = 1 outside the coUision manifold, 
i.e. on the domain r > 0. The energy surface Tih^ is of a shape of a hemisphere. The 
orbits are given by the intersections of Ti^ q with the angular momentum surfaces 
(not represented here). Easily noticed, in this case all orbits are periodic. 

If < a < 2, (yf(r) is positive for r < r* and negative for r > r* where 



Consequently, h{r) increases for r < r* and decreases for r > r*, and since 
h{0) > 0, the equation h{r) = has an unique solution. If a > 2, g{r) < 
thus h'{r) < 0, h{r) is a monotone decreasing function and h{r) = has a 
unique positive solution. 

Hence, for any a > the reduced smoothed flow admits two relative equi- 
libria located at Rf = (r^, 0, ±u^). 

3.4 Symmetries 

Observe that the plane u = is an invariant manifold. The dynamics is 
symmetric with respect to transformations of the form: 



The invariance under this symmetry implies that if 7(r) is a solution of (5) 
then also S'(7(r)) is a solution. We have the following: 

Proposition 3.3 If an orbit crosses the plane v = in one point it is S— symmetric. 
If it crosses the plane v = in two (distinct) points (with r > 0) it is 
S— symmetric and periodic. 

Proof The flrst part of the theorem follows from the uniqueness of solutions. 
The second part follows since the solution is 5*— symmetric and closed (and 
does not intersect a critical point). □ 




S : (r,f,M,r) {r,-v,u,-T). 
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Recall that for each fixed level h < of the energy, the reduced flow is 
constrained to the two dimensional surface E/j ^- Since the angular momen- 
tum is conserved, the orbits are determined by the intersection of S/j^^ with 
the two-dimensional angular momentum surfaces (see Figure 1). In other 
words, having e > fixed, to each pair energy-angular momentum pair {h, c) 
it corresponds a unique orbit 

7^^:=E;,,enFe. 

Using Proposition 3.3, it is sufficient to study the intersection of each orbit 
with the plane v = 0. Note that, on a given orbit, points of the form (r, 0,-u) 
with r > are turning points. Thus, for a periodic orbit, the two cuts with 
V — are the points where the relative distance r in between particle at- 
tains its maximum and its minimum. For a non-periodic orbit, the cut v — 
corresponds to the maximal relative distance. 

Since the graph of E/^^^ is symmetric with respect to the horizontal axis 
ti = 0, we can restrict our domain to m > 0. Note that from (13) we have 
that M > if and only if c > 0. So, fixing h < 0, the orbits are given by the 
intersection of 

with the curves Uc{r) as defined by 

2 — a 

Uc{r)r^2~ = c if a < 2, c > 0, 

a — 2 

Ucif^) = cr^~ if a > 2, c > 0. 

If this intersection is void, then there is no orbit, whereas if the w/i,e(r) and 
Uc{r) intersect in one (non-tangential) point, then the corresponding orbit is a 
fall to/escape from collision. Further, if the Uh^^{r) and U(.{r) have two distinct 
points in common, then the corresponding orbit is periodic. Finally, if Uh^^{r) 
and Uc{r) are tangent to each other, then the point of tangency corresponds 
to a relative equilibrium. 

3. 5 Topological equivalence of the reduced flows 

Let e > and consider the intersections of M/i,o(r) and Uh,e{T^) with Uc{r), 
where c > 0. 

Case a < 2. For c > 0, all orbits are periodic, exception being the fixed points 
Rq and Rf. For c — 0, all the orbits leading to/ejecting from collision (see 
Figure 2). Note that outside the collision set c — 0, the real and smoothed 
fiows are topologically equivalent. 

Case a = 2. The angular momentum curves become horizontal lines Uc{r) = c 
(see Figure 3). All orbits of the real flow are either leading to/escaping from 
collision, whereas all orbits of the smoothed flow are periodic. The real and 
smoothed flows are not topologically equivalent. 
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Fig. 2. Intersections of Uhfi{r) (top solid curve) and Uh^e{f) (bottom solid curve) 
with Uc{r) (dashed curves). Case a <2. 




r 



Fig. 3. Intersections of Uhfi{r) (top solid curve) and Uh,e{f) (bottom solid curve) 
with Uc{r) (dashed curves). Case a = 2. 

Case a > 2. First, note that all orbits of the real flow are all leading 
to/escaping from collision (see Figure 4). 

To determine the nature of the smoothed flow, we have to solve Uh,e{T^) — 
Uc{r). That is J fe{f) — cr^ , and further 



If c = 0, the previous relation becomes 



(r2 + e2)"/2 

The positive solution is unique and it is given by := \J {l/\h\y/" — e'^. This 
solution corresponds to an ejection-collision orbit with no spin. 



= 0. 



(14) 
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Fig. 4. Intersections of Uhfi{r) (top solid curve) and Uh^e{f) (bottom solid curve) 
with Uc{r) (dashed curves). Case a>2. (a = 3). Mo 



If c 7^ 0, let us consider xij') '■= ^ 2 ^2 °' ■ "^^^ positive solutions of 

2 



x{r) = are given by the positive solutions of ^(r) := — i^-^j + + 
= 0. Note that ip{0) = e^, and ip{r) ^ cxo as r — > cxd. Also, ip{r) has a 

unique minimum attained for r = r^^ := I - ( ^ j 1 . Moreover 'i/j{rm) — 

4 

(r, / „ \2/a\ 2a— 4 / „ \2/o TO 
j (^) [2/a — l] + e''<0 provided e is sufficiently small. It 

follows that ip{r) has two zeroes and thus x(r) has two zeroes. 

Thus x(r) — 2h has two solutions when h < and h is larger than the 
minimum of x(r), one solution when h is equal to the minimum of x(r) and 
zero solutions otherwise. We deduce that for h in the domain of interest (that 
is for non-void intersections of Uh,e{i^) with Uc{r)) and for c 7^ 0, the orbits 
for the smoothed flow are periodic, provided that e is small enough. Since the 
real flow consists of only collision/ejection orbits, the real and the smoothed 
flows are not topologically equivalent. 



4 The flow of the smoothed amended potential 

In the previous section we have shown that for a > 2 the flow associated 
to the smoothed potential is not topologically equivalent to the real flow. To 
understand why this is the case we return to the initial set-up of the problem. 

Recall the Hamiltonian of the real flow 

Ho{x,y,Px,Py) = I (pI+pI) - / (15) 
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or in polar coordinates: 



Ho{r,9,Pr,Pe) = 




It is clear from the expression above that the behaviour near coUision is dom- 
inated by the term i/7-™-'J^{2,a}_ follows that for a < 2, modifying the po- 
tential only does not change the main features of the flow near collision. This 
is not the case for a > 2. Here, smoothing of the dominant term allows the 
centrifugal term Pg/r^ to dominate near collision and to change the character 
of the motion. 



Recall that for the motion of a point in a central field on the plane, the 
distance from the centre of the field varies as in the one dimensional problem 
with a potential : 



C2 



where c := pg{t) = const, is the conserved angular momentum (see, for in- 
stance [3]). The function V{r) is usually called the amended or effective po- 
tential. The regions of motion together with the orbits' type (i.e. bounded or 
unbounded) are determined by the inequality V{r) < h. 

When smoothing is applied, the preservation, as much as possible, of the 
character of the orbits is necessary. For the homogeneous interaction with 
a >2, this may be achieved by smoothing the amended potential: 



K(r) := 



2(r2 + e2) (r2 + e2)"/2' 



As readily seen from Figure 5, the orbits are of similar type in both non- 
smoothed and smoothed problems, as long as c is chosen so that K(0) < 0. 
The regions of motion are determined by K(0) < h < and so, to have a 
non-void orbit set, c must be chosen such that: 



(17) 




To decide the topological character or the orbits, we return to initial set-up 
and consider the smoothed Hamiltonian 

He{r,9,pr,pe) := ^ {pI + p^^) " + ,2)a/2 ^ «>2- (18) 

Now we perform a study similar to one in the previous section. Omitting the 
details, the main steps are: 
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Fig. 5. The non-smoothed amended potential Vq and the smoothed amended Ve 
potential for a > 2. 

- first, we apply transformations similar to (3) and (4): 



V = Pr{r + e j 4 ; 
and introduce a new time parametrization via: 

dr = r^~^ dt] 



(19) 



(20) 



the new vector field is regularized and is given by: 



a+2 



vr 2 (r^ + e 



T 2 [r 



2 (U^ 



. a+2 , 2 , 

U = ur 4 (r^ + e 



U' 



— 2\ " + 2 / 9 , On "+4 

. 2 )^ 4 (r^ + e^j 4 uv] 



(21) 



- since ^ does not appear explicitly in the {r,v,u) equations above, we reduce 
the phase space to three dimensions by factoring out the fiow hy S^; 

- in the reduced space (r, v, u) the energy surfaces take the form 

E,,, = {{r,v,u) I u' + v' = 2(l- \h\{r' + 

they are surfaces of revolution and, as it may be easily noticed, their topology 
is identical; 



the angular momentum surfaces are given by 
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(a) (b) 

Fig. 6. (a) Intersections of Uhfi{r) (top solid curve) and Uh,e{f) (bottom solid curve) 
with Uc{r) (dashed curves). Case a = 2. (b)Intersections of Uhfi{r) (top solid curve) 
and Uh,e{f) (bottom solid curve) with Uc{r) (dashed curves). Case a> 2. 

{{r,v,u)\u^ c(r^ + e^) }; 

-by symmetry, each orbits may be represented by its intersection with the 
plane v = 0. Therefore, it is sufficient to study the intersections of the curves 

unAr) = ^/W) ■■= V2(l-|/i|(r2 + 62)"/2), 

with 

Wc.eW := c(r^ + e^)^, c > 0. 

For « > 2 all orbits are non-periodic, leading to or ejecting from the collision 
manifold M, := {{r,v,u) | r = 0, + ^2 ^ 2 (l - |/i|(e2)"/2j}. For a = 2, the 
angular momentum curves Wc,e(^) are horizontal lines (see Figure 5), whereas 
for a > 2, they are increasing and have initial value iic,e(0) = ce^^ (see 
Figure 6). The admissible values for c are obtained by requiring iic,e(0) < 

a — 2 I 

w/j^g(0). Thus we obtain ce^~ < \/2(l — |/i|e"), which, after some algebra, 
becomes condition (17). More importantly, the real and smoothed flows are 
topologically equivalent. 

5 Conclusions 

In this paper we discussed the topological equivalence of the non-smoothed 
and the smoothed flows given by motion with a homogeneous potential of the 
form — 1/r", a > 0. The analysis was performed outside the coUision set. We 
deduced that for a < 2 the two flows are topologically equivalent, and showed 
that this is not the case for a > 2. For the latter situation, we introduced 
the idea of smoothing the amended potential and showed that, with such a 
modification, the two flows are topologically equivalent. 
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Smoothing of the amended potential in a two degree of freedom system 
might be considered a first approach to a more general problem: given a me- 
chanical system with symmetry and with a singular potential, what is the 
best way to apply smoothing? Our investigation shows that the presence of 
non-inertial terms has to be treated carefully. Significant modifications of the 
global orbital picture might appear due to interplay of the potential and cen- 
trifugal forces, especially when the system passes close to a degenerate (e.g. 
for A^-body problems, a collinear) configuration. These issues will be discussed 
elsewhere. 
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